Problems in suppressing cooling flows in clusters of galaxies by 

global heat conduction 



Noam Soker 



c ^-. ABSTRACT 

O 
O 

CN ■ I use a simple analytical model to show that simple heat conduction models 

i-^ ■ cannot significantly suppress cluster cooling flows. I build a static medium where 

■ heat conduction globally balances radiative cooling, and then perturb it. I show 

that a perturbation extending over a large fraction of the cooling flow region 
and with an amplitude of ~ 10%, will grow to the non-linear regime within a 
Hubble time. Such perturbations are reasonable in clusters which frequently ex- 
perience mergers and/or AGN activity. This result strengthens previous findings 
which show that a steady solution does not exist for a constant heat conduction 
(N| ■ coefficient. 
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J^j ■ 1. Introduction 

Recent Chandra and XMM-Newton observations resolved the inner structure of the 
cooling flow region in clusters of galaxies. Traditionally, this region is defined as the region 
where radiative cooling time is shorter than the age of the cluster; not all clusters possess 
cooling flows (e.g., Fabian 1994). While hot, T~3-8x 10 7 K, X-ray emitting gas seems 
to cool at rates of 10 — 500 M Q yr -1 , there are almost no indications for gas cooling below a 
temperature of ~ 10 7 K (e.g., Kaastra et al. 2001; Peterson et al. 2001, 2002; Fabian 2002). 
This renewed interest in an old idea (e.g., Bregman & David 1988), that heat conduction 
from the intra-cluster medium (ICM) outside the cooling flow region supplies the energy 
carried by radiation (e.g., Narayan & Medvedev 2001; Voigt et al. 2002; Fabian, Voigt, & 
Morris 2002; Ruszkowski & Begelman 2002). Bregman & David (1988) show that the heat 
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conduction needs a fine tuning in order to explain the temperature and density profiles in 
cooling flows, and that a solution with a constant heat conduction coefficient is unstable, 
i.e., no steady solution exist. Norman & Meiksin (1996) try to increase the heat conduction 
efficiency by allowing the cooling flow to stretch magnetic flux loops, which then reconnect; 
cooling flow is not completely suppressed, but the cooling rate is substantially reduced in 
their model. Other works argue, from different reasoning, for a heat conduction coefficient 
k much bellow the Spitzer value for unmagnetized plasma, such that heat conduction does 
not play a global role in clusters and cannot prevent cooling flows (e.g., Pistinner & Shaviv 
1996; Markevitch, Vikhlinin, & Forman 2002; Loeb 2002). 

In a recent paper Zakamska & Narayen (2003; hereafter ZN03) fit five cooling flow 
clusters with a simple hydrostatic model where heat conduction balances radiative cooling. 
The thermal conductivity is ~ 30% the Spitzer conductivity. They argue that this model is 
stable against local transverse perturbations. For five cooling flow clusters they could not 
fit such a model, as they require too high conductivity. In this paper I show that such a 
solution is unstable to large large-scale radial perturbation with amplitudes > 10% in the 
temperature, over the age of a cluster. Hence, either the heat conduction does not play 
a role, or the cooling flow is much younger than the cluster, as proposed in the moderate 
cooling flow model (Soker et al. 2001). In either cases, heat conduction does not help in 
explaining the large scale behavior in cooling flow clusters, although it may be important 
locally, i.e., conduction to cold clouds in the inner region. In §2 I build a simple hydrostatic 
model with a power-law ICM. In §3 I use this model to argue that the solution is not stable 
over a long time, hence strengthening the results of Bregman & David (1988). In light 
of the renewed interest in heat conduction, (sometimes ignoring previous results), I aim at 
reaching this conclusion via a simple model. My short summary is in §4. 



2. A Simple Steady State Hydrostatic Solution 

To show that a solution where conduction balances radiative cooling is unstable, I start 
by building a simple spherically symmetric model with a power-law medium 

P = P or -V; T = T rP; n = n r~ v -^ (1) 

where P, T, and n are the thermal pressure, temperature, and electron density, respectively, 
and r is a dimensionless radial coordinate. The profiles used in equation (1) imply that 
the treatment is correct for r > 0, or more practically for r > 0.1. This does not change the 
conclusions of the present paper, because the instability found here and most of the radiative 
energy loss occur at much larger distances from the center. Like ZN03 I take the cooling 



- 3- 



rate per unit volume and the heat flux to be j = Kin 2 T 1/2 , and F = -fK 2 T 5/2 dT/dr, 
respectively, where K 1 and K 2 are constants, and / is the ratio of the conductivity to the 
Spitzer value. The total radiated power, a free-free emission, inside a radius r is given by 

£rads = f K x n 2 T x l 2 4ixx 2 dx = M K^Tl' 2 (<o - 477 - 3/5)~ 1 r 3 - 2 "-i /3 , (2) 
Jo 

where subscript 'S' stands for the steady state solution. Most of the energy is radiated just 
inward to r, such that the detailed structure close to the center is not important. The rate 
of inward heat flow across a spherical surface at r is 

Cheats = -4vrr 2 F = A-KfK 2 (3Tl /2 r^ +1 . (3) 

In a steady state hydrostatic model Cheats — -^rads for each radius, which implies two equali- 
ties: the equality of the powers of r in equations (2) and (3) and the equality of the constant 
terms in these two equations. The first equality gives 

ri = 1 - 2.5/?, (4) 

and the second equality gives 

fK 2 f3T* = 2K in 2 (6 -4 V - 3/3)- 1 . (5) 

From the solutions of ZN03, I find rj ~ 0.3. Equation (5) implies that when / is changed 
by a factor a, this may result in a modest change in the temperature and a change in the 
density by a factor of ~ a 1 / 2 ; this is compatible with the finding of ZN03. 

3. Perturbing the Temperature Profile 

I take a simple form for the large scale temperature radial perturbation, AT — Ts{r — \), 
such that the new temperature profile is 

T / = T + AT = T r /3 + T (5 (r-l), (6) 

where T is given in equation (1), and the derivative of the temperature is (for r > 0) 

^ = pry-' + t s . (7) 

For simplicity I assume that the pressure profile does not change, hence the new density 
profile is 

ri = n(T/V). (8) 
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A full accurate treatment should consider the perturbation in pressure as well. The pressure 
perturbation is neglected here for the following reasons. (1) As stated in Section 1, the aim 
of the present paper is to present a simple treatment in order to clearly demonstrate the 
global radiative instability. (2) This assumption is accurate for a case of constant pressure. 
Although the pressure in clusters is not constant, it is still quite shallow. As derived below, 
the instability which occurs for a constant pressure is not much different from the one for a 
shallow pressure profile. (3) Although being simple, the result of the present treatment is in 
accord with the numerical results of Bregman & David (1988). This makes the treatment 
trustworthy. 

Only the first order in T$ is retained. For example, 



(T') 5/2 = [Tor' 3 + T 5 (r - l)f 2 = T 5/ V^ 2 



1 + 



r-l) 



+ 0(T|). 



(9) 



Keeping only the first order is not justified very close to the center where the unperturbed 
temperature is very low. As noted earlier, this does not affect the conclusions here, because 
the instability occurs farther out in the cluster, and because in real clusters the very central 
region is expected to be dominated by the AGN activity any how. Substituting equations 
(7) and (9) in the expression for the total rate of inward heat flow across a spherical surface 
of radius r, similar to equation (3) for the steady state heat flow LheatSj gives to first order 
inT 5 

T s (I 5 5" ~ 

- I — r A r 

2 
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Expanding (n') 2 (T') l l 2 to first order in Ts, and integrating, similar to equation (2) for the 
steady state cooling rate Z/ ra dS, and using equation (4) for rj, gives 



7 rad 



>radS 



2 + 7/? _ 2 + 7/3 

2 T r? \2 + 5P ~ 4 + 5/? ? 



(ii) 



The net cooling rate of the region inward to r is given by L rad — L hcat , which result in 
a cooling time r coo i(r). Let r rad (r) be the cooling time calculated when only the radiative 
cooling L ra ds is considered, i.e., the "standard" cooling time at radius r. (The cooling time 
calculated at radius r, is not exactly the same as the average cooling time of the gas inward 
to r. However, it was noticed above that most of the radiative cooling rate inward to r 
occurs close to r. I therefore neglect the difference between these two values.) The ratio 
between r ra( j(r) and r cool (r) is given by 



r rad( r ) _ -^rad ~ Aicat 
Tcool (r ) -^radS 



(12) 
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Substituting for Lhcat and L ra d from equations (10) and (11), respectively, and using Lhcats = 
-^radS) gives to first order in T$ 

r rad (r) T s {6 + 21P 5 6 + 21(3 1 5 \ 



rcooi(r) T(r) v ' T(r) V r-1 

TUr) - T ' -(4.3 - 6.2r) = ^ f 4 - 3 - 6 - 2r ) , for = 0.4. (15) 



T-cooi(r) T(r)V4 + 10/? 2 8 + 10/3 /3 2 
As mentioned above (5 ~ 0.3. Equation (13) gives then 

and 

Tcooi(^) T(r) v "~ ~' ^(r) V r ~l 
Note that f3 = 0.4 is the case of constant pressure. Although this is not realistic, in the 
assumption made here of neglecting the pressure variation, the treatment of this case, where 
there is no pressure gradient, is accurate. The last two equations show that there are no 
fundamental differences between the instabilities in the case of constant pressure and the 
case of a shallow pressure gradient appropriate for clusters. 

Let us examine the nature of the instability. In a case of no radiative cooling, or a 
slow radiative cooling, the restoring mechansim of a temperature perturbation is the heat 
conduction itself. Namely, when a region gets cooler the temperature gradient gets steeper, 
the energy flow via heat conduction from warmer to cooler regions increases, restoring the 
cooling region back to the equilibrium temperature. What was showen above is that under 
the conditions which exist in cooling flow clusters, when a large inner regions gets cooler, 
hence denser, the energy lose rate via radiative cooling, i.e., the radiative power, increases 
more than the heating power via the heat conduction does. For local perturbations, the 
short scale means steep temperature gradient, hence efficient heat conduction. This is why 
ZN03 find their perturbation to be stable. I instead consider large scale perturbations, where 
radiative energy lose occurs in a large region, as appropriate for the inner cooling flow region. 
Such large perturbations are reasonable in the violent environment near the cD galaxy in 
cooling flow clusters. The time evolution of the departure from the initial perturbative state 
occurs on the cooling time scale or longer (since some heat conduction may exist). Therefore, 
it is important (during the cluster age) only in the inner regions of cooling flow clusters (and 
only in cooling flow cluster) where cooling time is very short. 

For an instability to develop, the expression inside the parenthesis in the far right-hand 
side of the last two equations should be negative. This occurs for r > 1. However, for a large 
scale perturbation this occurs at relatively large radii, where the standard cooling time r ra d 
is long. There is also an instability at much smaller radii. For these, the perturbation, as 
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chosen here, must be a large scale perturbation: It starts at r = 1 (eq. 6), and the instability 
occurs only for r < n = 4.2/8.52 = 0.5 for (3 = 0.2, and r < n = 4.3/6.2 = 0.7 for (3 = 0.4. 
I consider the case (3 = 0.2 (the results are not sensitive to the value of /3), at r = 0.25. For 
example, the perturbation starts at r = 1, which in real units I take at R = 60 kpc, and 
I consider the region around r ~ 0.25, or in real units R ~ 10 — 20 kpc. The "standard" 
cooling time in this region in cooling flow clusters is r ra d(i? = 10 — 20 kpc) ~ 2 — 5 x 10 8 yr. 
Hence the cooling time, scaled with a perturbation of ~ 10% is 

for a perturbation given by equation (6) with r = 1 at R ~ 40—80 kpc. For these parameters, 
the perturbation will grow to the nonlinear regime during the age of the cluster, ~ 10 10 yr. 
Therefore, a model with no cooling flow requires, either that (1) the cooling flow is much 
younger than the cluster, e.g., there is an episodic heating of the cooling flow region on a 
time scale of ~ 10 9 yr (e.g., Soker 2001), or (2) there is another continuous energy source, 
e.g., an AGN at the center of the central cD galaxy (Ruszkowski & Begelman 2002). In 
either cases, heat conduction by itself cannot "kill" the cooling flow. 

The following should be noted. The instability occurs for r < r iy whereas the region 
r\ < r < 1 is stable. However, this stable region cannot stabilize the unstable region r < r\. 
This is because the instability evolves on a cooling time scale, which is much shorter in the 
inner unstable region than in the stable region. The non-existence of a steady state solution 
was found by Bregman & David (1988) in their numerical calculations. 



4. Summary 

New results by the Chandra and XMM-Newton X-ray telescopes, which indicate a very 
low mass cooling rate at temperatures bellow ~ 1 kev (see review by Fabian 2002), renewed 
interest in heat conduction from the outer ICM as an energy source to compensate for ra- 
diative cooling at lower temperatures. The goal of the present paper is to use a simple 
analytical model to show that (simple) heat conduction models cannot significantly suppress 
cluster cooling flows. I used a power-law ICM and built a static medium where heat conduc- 
tion globally balances radiative cooling, and then perturbed it. I showed that a large scale 
temperature perturbation with an amplitude of ~ 10%, will grow to the non-linear regime 
within a Hubble time. The perturbation should extend over a large fraction of the cool- 
ing flow region. Such perturbations are reasonable in clusters which frequently experience 
mergers and/or AGN activity, e.g., X-ray deficient bubbles (McNamara 2002). This result 
strengthens the finding of Bregman & David (1988), who use numerical calculations to show 
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that a steady solution does not exist for a constant heat conduction coefficient. 

This instability over a time scale of ~ 1 9 yr, implies that another energy source should 
either make the cooling flow a short-lived episodic process, or that it continuously supplies 
energy. In either cases the heat conduction does not seem to play a significant role; the 
major role is played by the other energy source(s). A popular energy source is AGN activity, 
which was proposed to reduce the mass cooling rate in cooling flows in galaxies (e.g., Binney 
& Tabor 1995; Ciotti & Ostriker 1997, 2001; Jones et al. 2002), and in clusters of galaxies 
(e.g., Soker et al. 2001; Churazov et al. 2002) 

To the instability found in the present paper, I add the failure of ZN03 to build a steady 
model with heat conduction to five clusters, the claim made by Sun et al. (2003) that 
generally heat conduction can't compensate for radiative cooling in cluster cooling flows, 
the indications that heat conduction is suppressed at other regions of clusters (Loeb 2002; 
Markevitch et al. 2002; Nath 2003), and the required fine tuning argued by Bregman & 
David (1988), and conclude that heat conduction cannot explain the properties of cooling 
flow clusters. 

I thank the referee and Larry David for useful comments. This research was supported 
in part by grants from the US-Israel Binational Science Foundation and the Israel Science 
Foundation. 
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